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EXCHANGEABLE, STATIONARY AND ENTANGLED CHAINS OF 

GAUSSIAN STATES 

K. R. PARTHASARATHY AND RITABRATA SENGUPTA 


Abstract. We explore conditions on the covariance matrices of a consistent chain of mean 
zero finite mode Gaussian states in order that the chain may be exchangeable or stationary. For 
an exchangeable chain our conditions are necessary and sufficient. Every stationary Gaussian 
chain admits an asymptotic entropy rate. Whereas an exchangeable chain admits a simple 
expression for its entropy rate, in our examples of stationary chains the same admits an integral 
formula based on the asymptotic eigenvalue distribution for Toeplitz matrices. An example of 
a stationary entangled Gaussian chain is given. 
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1. Introduction 

The importance of finite mode Gaussian states and their covariance matrices in general quan¬ 
tum theory as well as quantum information has been highlighted extensively in the literature. 
A comprehensive survey of Gaussian states and their properties can be found in the book of 
Holevo [Holll]. For their applications to quantum information theory the reader is referred to 
the survey article by Weedbrook et al [WPGP’^12] as well as Holevo’s book [Holl2]. For our 
reference we use [ADMS95, ParlO, Paris] for Gaussian states and for notations in the following 
sections we use [PS15]. While working on this paper, Mozrzymas, Rutkowski, and Studzihski 
had posted in arXiv [MRS15] an article similar in spirit, but for finite dimensional Hilbert 
spaces. 

In the present paper our concern is with a chain of finite mode Gaussian states constituting a 
consistent sequence exhibiting properties like exchangeability, stationarity, and entanglement. 
All these notions are easily translated into properties of infinite covariance matrices. 

If p is a state of a quantum system and Xj, z = 1,2 are two real valued observables, or 
equivalently, self-adjoint operators with finite second moments in a state p then the covariance 
between Xi and X 2 in the state p is the scalar quantity Tr p—TrXip ■ TrX 2 p, which 

is denoted by GoVp(Xi,X 2 ). Suppose q'i,pi; q 2 ,P 2 ] ■ ■ ] QmPn are the position - momentum 
pairs of observables of a quantum system with n degrees of freedom obeying the canonical 
commutation relations. Then we express 


(Xi,X 2 ,-- - ,X 2 n) = {qi,Pl,q2,P2, - ■ ■ ,qn,Pn)- 
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If p is a state in which all the Xj^s have hnite second moments we write 
(1.1) [[Cov,(X„X,.)]], 1,J e {1,2,-- - ,2n}. 

We call Sp the covariance matrix of the position momentum observables. If we write 


( 1 . 2 ) 


Jin — 


0 1 

-1 0 


0 1 

-1 0 


0 1 

-1 0 


or equivalently 0” 



1 

0 


for the 2n x 2n block diagonal matrix, the complete Heisenberg 


uncertainty principle for all the position and momentum observables assumes the form of the 
following matrix inequality 


(1.3) Sp + 2'^2n ^ 0. 

Conversely, if S is any real 2n x 2n symmetric matrix obeying the inequality Sp + ^Jin > 
0, then there exists a state p such that S is the covariance matrix Sp of the observables 
5'DPij Q 2 ,P 2 ] ■■■ ; QntPn- Ip such a case p can be chosen to be a Gaussian state with mean 
zero. In view of this we make a formal dehnition. 


Definition 1.1. A 2n x 2n real symmetric positive matrix S is said to be a G-matrix if it 
satisfies the inequality 

(1.4) S + 2'^2n > 0. 

Suppose 

(1.5) ^ = [[Afi], qje{l,2,---} 

is an inhnite matrix where each Aij is a 2A; x 2k real matrix and AJ = Aji for all i, j. For any 

hnite subset / = < f 2 < • • • < in} C. {1, 2, • • ■ }, let 

(1-6) S(J) = [KJ], r,5e{l,2,---} 

be the 2kn x 2kn matrix obtained from S by restriction to its rows and columns numbered 

ii < ii <■■■< in- 

Definition 1.2. We say that T, is a G-chain of order k if T,{I) is the covariance matrix of a 

kn-mode zero mean Gaussian state p{I) in the boson Fock space r(C^"') = ® 

where FLj denotes the j-th copy of the Hilbert space FL = r(C^), j = 1, 2, • ■ ■ . 

If / = {H < f 2 < • ■ • < in} and /' = {H < *2 < ■ ■ • < in+m} then clearly p(/) is the I- 
marginal of the state p{F)- Thus S describes a consistent family of zero mean Gaussian states 
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with finite number of modes which are multiples of k. It is a standard result of Gaussian states 
that E is a G-chain of order k if and only if the following matrix inequalities hold: 

(1.7) S({1,2, • • • ,n}) + > 0, n = l,2, 


Definition 1.3. IFe say 

and there exist two 2k x 

that S in equation (1.5) is an exchangeable G-chain if it is a G-chain 
2k matrices A, B such that 

(1.8) 

Aij — 1 

f H if j >i, 

,71 if j = i, 

[ if J < i. 

In such a case we write 

E 

= ^{A,B), 

(1.9) 

E(/) 

= W;A,b). 


It is clear that for any two finite subsets I and /' of {1, 2, • • ■ } with the same cardinality an 
exchangeable G-chain S satisfies 

and therefore the corresponding Gaussian states p(/) and p{I') are same, i.e. the quantum 
version of de Finitti exchangeability property [de 29] holds. 

In this paper we shall show that two 2k x 2k matrices A and B determine a G-chain E(y4, B) 
if and only ii B = 5 > 0 and A — B is the covariance matrix of a fc-mode Gaussian state. 

For any finite mode Gaussian state p with covariance matrix C denote its von Neumann 
entropy S{p) by S{C). we shall prove that for any exchangeable G-chain E(y4, B), the sequence 
{7S'(E({1, 2, • • ■ ,n},74, 5))} decreases monotonically to the limit S{A — B). 

Definition 1.4. We say that a G-chain S given by (1-5) is stationary if there exist 2k x 2k 
matrices A, Hi, H 2 , • • • such that 

( ^ 

(1.10) Ai,- = { Bj_i 

I 

for all i, j. 

For such a stationary G-chain the Gaussian states p{I) and p{I 1) are same for any finite 
subset J C {1, 2, • ■ ■ }. This translation invariance property shows that our definition of station- 
arity is similar to such a notion in classical theory of the stochastic processes [Doo53]. It is an 
interesting problem to find necessary and sufficient conditions on the matrices A, Hi, H 2 , • • • so 
that (1.10) yields a stationary G-chain. Even though we do not have an answer to this question 
we shall construct a large class of examples of such stationary G-chains. 

For any stationary G-chain S, consider the von Neumann entropy of the Gaussian state 
p({l, 2, • • • , n}. A, H), namely, S'(S({1, 2, • • ■ , n} : A, H)). The strong sub-additivity property 
of entropy [LR73] implies that the sequence {^^(Edl, 2, ■ ■ ■ , n} : A, H))} decreases monoton¬ 
ically to a limit *F(S), which may be called the entropy rate of the stationary chain. Using 
the Kac-Murdock-Szego theory for asymptotic eigenvalue distributions of Toeplitz matrices 


if i = J, 
if j > i, 
if j <i 
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[KMS53, Tro84] we shall get an integral formnla for the entropy rate of an interesting example 
of a stationary G-chain of order k. 

When A; = 1, we shall constrnct a stationary G-chain in which states of the form p({l,2}) 
and hence p({l, 2, • • • , n}), n = 2, 3,4, • • • are entangled. However it wonld be more interesting 
to find examples of stationary G-chains in which, for some fixed finite nonempty snbset I of 
{1, 2, • • • }, states of the form p{I U (/-|-n)) with marginals p{I) remain entangled for arbitrarily 
large n. That wonld ensnre preservation of entanglement after an arbitrarily large lapse of time. 

2. Exchangeable G-chain 

Onr hrst main resnlt gives necessary and snfficient conditions on a pair [A, B) of real 2k x 2k 
matrices so that the inhnite matrix T,{A,B) in Definition 1.3 is an exchangeable G-chain. 

Theorem 2.1. Let {A,B) be a pair of real 2k x 2k matrices. Then I1{A,B) is a G-chain if 
and only if A and B are nonnegative definite and A — B is a G-matrix. 

Proof. Ghoose and £x a positive integer n and pnt I = {1,2,-- - ,n}. Write 

Denote by the identity matrix of order n and the npper triangnlar nilpotent matrix with 
all its npper triangnlar entries eqnal to 1 and the rest eqnal to zero. Let 

= .ly. 


a nnit colnmn vector of length n. Then 

In + Nn + N'^ = n |V’n)(V^n| • 

Then we have 

B) = A®In + B®Nn + B'^ ®Nl. 
Noting that J 2 kn = J 2 k ® In we can now write 


I^n^A, B) + ^J2kn — 


A -h -J2k - -(-B + 


+ ( ^ + 


5 {In - \ fjn){i’n\) 
1){B + ) 0 


( 2 . 1 ) 


+ -{B-B^)®{Nn-N^). 


In order that T{A,B) may be a G-chain it is necessary and snfficient that the left hand side 
of (2.1) is nonnegative dehnite for every n. To prove necessity we mnltiply hrst both sides of 

(2.1) by In 0 \‘ipn){'f’n\ and take relative trace over the second component. Noting that for any 
real vector If)) in the second Hilbert space, (i/'lW — Nnl'^) = 0 we get the ineqnality 

(2.2) A —J 2 k — 1)(-B T -B”^) ^ 0 Vn. 


Dividing by n — 1 and letting n —)■ cx) we get 


(2.3) 
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Choosing an arbitrary real unit vector \'ip) in the range of — |'0„)('0„|, multiplying both sides 
of (2.1) by / ® |'0)('0| and tracing out over the second Hilbert space we get 

(2.4) yl + lj2fc-^(5 + H^) >0. 

In other words A — \{B + B^) is a G-matrix. Now we consider the complex unit vector 

I0n) = ^[1,0;,- 


where a; = e " , an n-th root of unity. Simple algebra shows that 

1 


{(t)n\Nn\(t)n) = 


and therefore 


{(pnWn - iVn I0n) = 2*/m 


U — 1 
1 


bJ 


TT 

- = I cot —. 
1 n 


Now, multiplying both sides of (2.1) by ® |0n)(0„| and tracing over the second Hilbert space 
we get the inequality 


(2.5) 


7 1 B — TT 

A + -J2k --(B + B^)+ i- -cot - > 0 

2 2^ ’ 2 n- 


for n = 1, 2, • • ■. Multiplying by tan ^ (for n > 3) and letting n —)• cx) we get the inequality 

^-{B-B^)>0. 

Since the left hand side is a Hermitian matrix with trace zero it follows that B = B^ and (2.4) 
implies that A — B is a. G-matrix. Together with (2.3) the proof of necessity is complete. 

To prove sufficiency, observe that A — B being a G-matrix and B being positive implies that 
A — B + nB is a G-matrix for every n = 1,2, • • •. Now the identity (2.1) implies that 

'^n{A, B) + -J2kn = ~ B) -\- ^'^2fc 

for every n and therefore S„(y4,i?) is a G-matrix for every n. In other words E(y4,i?) is a 
G-chain. □ 


(4 - |'0n)('0n|) + A+ {n-l)B + -J2k ® \'tpn){'tpr. 


Corollary 2.1. In any exchangeable G-chain Yi{A,B) of order k, for every finite set I C 
{1, 2, • ■ ■ the underlying Gaussian state p{I) is separable. 

Proof. Without loss of generality we may assume that I = {1,2,-- - ,?t,} for some n. Then the 
covariance matrix of p{I) is equal to the n x n block matrix 


'A 

B •• 

• B' 


T 

1 ■■ 

■ 1' 

B 

A ■■ 

■ B 

= (A- B)(g)In + B^ 

1 

1 ■■ 

■ 1 

_B 

B 

■ A 


.1 

1 ■■ 

■ 1. 


By Theorem 2.1, (H — H) G) 4 is the covariance matrix of an n-fold product Gaussian state 
and the second summand on the right hand side of the equation above is a nonnegative definite 
matrix. Hence by Werner and Wolf’s theorem [WWOl] p{I) is separable. □ 
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3. Examples oe stationary G-chain 


Let A, B be real 2k x 2k symmetric matrices. For any fixed j = 1,2, ■ ■ ■, denote by Ai^{A, B) 
the infinite block matrix all of whose diagonal blocks are eqnal to A, {n,n + j)-th and {n+j, n)- 
th blocks are eqnal to B for every n and all the remaining blocks are zero matrices of order 
2k X 2k. For example, 


A\A,B) 


'ABO 0 0 
B A B 0 0 
0 B A B 0 


Denote by A^(y4, 5) the 2kn x 2kn matrix obtained by A^{A,B) by restriction to the first n 
row and colnmn blocks. For example, 


Al{A,B) 

Al{A,B) 


A B 
B A 

'A 0 B' 
0 A 0 
BOA 


and so on. 

Onr first resnlt gives a necessary and snfficient condition for A^{A, B) to be a G-chain. 


Theorem 3.1. Let A, B be a pair of 2k x 2k real symmetric matrices. In order that A^{A, B) 
may he a G-chain of order k it is necessary and sufficient that A-\-tB is a G-matrix for every 
f e [-2,2]. 


Proof. Denote by Lf the npper triangular matrix whose (j -|- l)-th upper diagonal entries are 
all equal to 1 and all the remaining entries are zero. Thus Lf is defined for 1 < j < n — 1. 
Then 


(3.1) Ai(A, B) = A»I„ + B» (Li + (Lif). 

Consider the spectral decomposition of the n x n symmetric matrix 

n 

(3.2) Li + {Lif = Y, An. \f^nr){i^nr\ 

r=l 

where {A„r : = 1, 2, ■ • • , n} are the eigenvalues and {iV'n.) : = 1, 2, ■ ■ ■ , n} are the cor¬ 
responding orthonormal basis of eigenvectors for Li -|- (Li)'^. Since each Li is a matrix with 
operator norm equal to unity and therefore Li -|- (Li)'^ has operator norm not exceeding 2 it is 
clear that 


(3.3) lAfir-l <2, 1 < r < n, n = 1, 2, • • • . 

Equations (3.1)~(3.2) imply 

n 

(3.4) B) + -J2kn = ^nrB 2 ® \'f’nr){'f’nr 

r=l 
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Thus A^(y4, B) is a G-matrix if and only ii A + XnrB is a G-matrix for each r = 1, 2, ■ ■ • , n. 
This together with (3.3) already proves the sufficiency part of the theorem. 

To prove necessity, we appeal to the theorem of Kac, Murdock and Szego [KMS53] . Gonsider 
the probability distribution 

1 

n ^ 

r=l 

where f ,?t, are as in (3.2). The left hand side of (3.2) is a Toeplitz matrix 

of order n for each n. Thus Kac, Murdock, Szego theorem implies that the sequence {/Xn} 
converges weakly as n —)■ cxd to the probability measure Lh~^ where L is the Lebesgue measure 
in the unit interval and h{t) = 2cos27rjt, t G [0,1]. This, in particular, implies that {A„r : 
r = 1, 2, • ■ ■ ,n, n = 1, 2, ■ ■ ■ } is dense in the interval [—2, 2]. The proof of necessity is now 
complete. □ 


Corollary 3.1. Let A, Bi, B 2 , ■ ■ ■ be real 2k x 2k symmetric matrices satisfying the condition 
that A + tBj is a G-matrix for every j = 1, 2, • • • and t G [—2, 2]. Suppose Pi,P 2 , ■■■, is a 
probability distribution on the set {1, 2, • • • }. Then the block Toeplitz matrix 


S(A;pi5i,P252,---) 


A piBi P 2 B 2 ■ ■ ■ 
PiBi A piBi P 2 B 2 
P 2 B 2 piBi A piBi 


is a stationary G-chain 

Proof. This is immediate from the theorem because 

00 

J:{A;piBi,P 2B2, ■■■) = '^PjA^{A,Bj) 

i=i 


and each A^{A,Bj) is a G-chain. 


□ 


4. Entropy rate of the stationary G-chain 

Suppose E = S(y4, Bi, B 2 , ■ ■ ■) is a stationary G-chain. For any G-matrix C denote by S{C) 
the von Neumann entropy of a Gaussian state p with covariance matrix C. Let 




S({1,2,--- ,n}), 

^(Sn). 


Proposition 4.1. The sequences {S'„ —monotonically decrease to the same limit 
S > 0 as n ^ 00 . Furthermore, Sn > Sn-i for all n. 


Proof. Gonsider three Gaussian quantum systems P, Q, R so that PQR is in the mean zero 
Gaussian state p({l,2, ,n-|-l}), Q in p({2, • • • ,n}), PQ in p({l,2, ,n}) and QR in 
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p({2, • • • , n + 1}). Then using stationarity we have 

S{p{PQR)) = Sn+l 

SipiPQ)) = Sn 

SipiQR)) = Sn 

S{pm = Sn-l. 

By the strong subadditivity property of entropy [LR73] we have 

Sn+l + Sn-1 < “^Sn 

or 

Sn+l ~ Sn ^ Sn — 

S incG 

_ {Sn - Sn-l) + {Sn-l - Sn- 2 ) + ' ' ' + (>^1 „ 

n n 

where Sq is dehned to be zero, it follows that ^ decreases monotonically to a limit ^ > 0. 
This also implies that Sn — Sn-i cannot decrease to —oo and hence Sn — Sn-i also decreases 
monotonically to S. This also shows that Sn > Sn-i for all n. □ 

Definition 4.1. ITe denote the limit S in Proposition 4S by *S(S) and call it the entropy rate 
of the stationary G-chain S. 

Theorem 4.1. Let S = he an exchangeable G-chain. Then 

S{T) = S{A-B). 

Proof If C, D are two G-matrices so is G © D and S{C ® D) = S{C) + S{D). By Theorem 2.1 
and equation (2.1) we have the identity 




Tn{A, B) = {A- B)® {In - |V’n)(V’n|) + {A P U - IB) ® \lpn){'<Pn\ 

where we have adopted the notations in the proof of Theorem 2.1. Thus 

(4.1) Sn = S{Tn{A, B)) = {n- l)^(Al - B) + S{A + fG^B). 

Denote by p^ the zero mean Gaussian state with covariance matrix A for any G-matrix A. 
Thus we have 


P 


A+n—lB 




where ^ ©^ 2 , bT(^) is Weyl or displacement operator at + i $,2 and 0(^) is the Gaussian 

density function with mean zero and covariance matrix (n — l)i?. By Proposition 6.2 of [OP93] 
we have 


S{APn-lB) = S{p^+^-^^) 


(4.2) 


< 


S{A)m^i + H{ct>) 


where H{(f)) is the Shannon differential entropy of the density function <f. Since (by [GT06]) 

(4.3) 


H{(f)) = k\og2'Ke + - logdet(n — IB) 
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it follows from (4.1)-(4.3) that 


n 


n 


< + — log 2716 + — login — 1)^^ det B 

n n 2n 


1 

< - 
n 


S (A) + k log 27re + ^ log det B 


k 

+ - login - 1). 
n 


Letting n —)■ cxd we get 


^(S) = S{A-B). 


□ 


Theorem 4.2. Let pi,P 2 , ■ ■ ■ be a probability distribution over {1,2, 3, • • ■ and let A and B 
be 2k X 2k symmetric real matrices satisfying the condition that A + tB is a G-matrix for every 
t G [—2,2], Let S be the stationary G-chain defined by the infinite block Toeplitz matrix 


S 


■ A piB P 2 B ■ ■ ■ 
PiB A piB P 2 B 
P 2 B piB A piB 


Then the entropy rate of S is given by 


S{E)= / S{A + h{s)B)ds 


where 


h{s) = 2'^^ pj cos 271 js, s G [0,1]. 

i=i 

Proof. We can express as 

= 21 ® 4 + 5 0 TM 

where is S restricted to its hrst n row and column blocks and Tn{p) is the Toeplitz matrix 
given by 



■ 0 

Pi 

P 2 ■ 

■■ Pn -1 


Pi 

0 

Pi ■ 

■■ Pn -2 

Tn{p) = 

P 2 

Pi 

0 ■ 

■ ■ Pn -3 


Pn-l 

Pn -2 

Pn -3 ■ 

0 . 


Let A„i, Xn 2 , • ■ ■ , Xnn he the eigenvalues of T„(p) and let |'0ni), \'f’n 2 ), •'' ) If’nn) the correspond¬ 
ing eigenvectors constituting an orthonormal basis for so that 

n 

+ XnjB) 0 \i>nj){i>nj\ ■ 

i=i 
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This shows that 


1 1 

-S{En) = -y^S{A + XnjB) 
n n 

j=i 

= j S{A + sB) d/i„(s), 

where fin is the probability measure defined by 

1 ” 

■ 

i=i 

By Kac-Murdock-Szego theorem fin converges weakly as n —)■ cx) to the distribution Lh~^ where 
L denotes the Lebesgue measure in [0,1] and 

OO 

h{s) = 2 E Pj cos 271 js. 

J=1 

Note that ||T„(p)|| < 2 and the eigenvalues Xnj lie in the interval [—2,2], Furthermore, the 
symplectic spectrum oi A + sB is a continuous function of s and hence the entropy S'(A + sB) 
is a continuous function of s in [—2, 2], Thus 

lim -5(S„) = / S{A +sB)Lh-^{ds) 

n—>-oo n J _2 

= f S{A + h{s)B)ds. 

Jo 

□ 

5. Entanglement in a stationary G-chain 

We have already noted in Corollary 2.1 that in any exchangeable G-chain T,{A, B) of order k 
all the underlying Gaussian states p{I) with > 2 are separable. In the class of examples of 
stationary G-chains in Corollary 3.1 it is natural to examine the presence of entanglement. It 
would be particularly interesting to construct an example of a stationary G-chain for which the 
underlying Gaussian state p(/U(/ +n)) with marginals p{I) and p{I + n) = p(/) is entangled for 
arbitrarily large n. That would imply the preservation of entanglement after an arbitrarily large 
lapse of time. For such problems we do not have any answer. However, we shall examine the 
special case of a stationary G-chain of order 1 where the matrices A, Bi, B 2 , ■ ■ ■ in Corollary 
3.1 are given by 


A 

Bj 


XI2 

B = b 


1 0 
0 -1 
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■A/2 

PiB 

P2B 


PiB 

XI2 

PiB 

P2B 

P2B 

PiB 

XI2 

PiB 


where A and b are positive scalars with A > |. We start with two elementary lemmas. Let 

■A/2 p] 
mB A 

(5.1) S = 

where pi,p 2 , - ■ ■ is a probability distribntion over {1,2, • • ■ } and A, /? as above. 

Lemma 5.1. The infinite block matrix S in (5.1) is a stationary G-chain of order one if 

Proof. By Corollary 3.1, S is a stationary G-chain of order one if AJ 2 + tB is a. G-matrix for all 
t G [—2, 2]. This is fnlfilled if the matrix ineqnalities 


\ + tb 

_ 

2 


> 0 . *€[- 2,21 


hold. Clearly this is satisfied if 6 < ^ (A^ ~ I) ^ • 
Lemma 5.2. Let A > |, c > 0. Then the matrix 


□ 


r = 


A 0 c 0 
0 A 0 -c 
c 0 A 0 
0 -c 0 A 


is the covariance matrix of an entangled 2-mode Gaussian state if 


2 V 4 


Proof. By Simon’s criterion [SimOO], T has the reqnired property if the following two matrix 
inequalities hold: 


(5.2) 

(5.3) 


<< 

1_ 

_c2 

'1 0 ■ 


1 

tOl'S* 

-1 

'1 o' 

L-i 


0 -1 


L-i 


0 1 


_c2 

'l 0 ■ 


[A -1] 

-1 

'l o' 

[-1 


0 -1 




0 1 


> 0 , 

> 0 . 


Simple algebra shows that (5.2) holds whenever < A^ — \. Inequality (5.3) reduces to 

\{l-d^)X |(l + d2)l 

[-1(1+ d^) {l-d^)X 

where 


0 
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We now choose d < 1 such that the determinant of the matrix on the left hand side of (5.4) is 
strictly negative. This dually leads to the inequality c > A — Thus the inequality A-i< 

c<^ — I is sufdcient to ensure entanglement. □ 


Proposition 5.1. Let | <A< A — | < b < \ - Suppose pjb > A — | for some j. 

Then the 2-mode Gaussian state determined by the stationary G-chain S defined by 

(5.1) is entangled. 


Proof. Since ^ < A < |, the interval (^A — | (A^ ~ ^ j is open and nonempty. Thus it is 

possible to choose b in this interval. By Lemma 5.1 the block matrix E in (5.1) is a stationary 
G-chain of order one. The covariance matrix of the 2-mode state p({l,j}) dedned by the 
G-chain E is equal to 

A 0 Pjb 0 

0 A 0 —Pjb 

Pjb 0 A 0 

0 —Pjb 0 A 

where by hypothesis 

By Lemma 5.2 it follows that /9({l,j}) is entangled. □ 


Remark 5.1. It follows from Proposition 5.1 that for ^ <A<|, A-|<6<^(A^ 
matrix 


■ A 

0 

b 

0 

0 

0 




0 

A 

0 

-b 

0 

0 




b 

0 

A 

0 

b 

0 

0 

0 


0 

-b 

0 

A 

0 

-b 

0 

0 


0 

0 

b 

0 

A 

0 

b 

0 


0 

0 

0 

-b 

0 

A 

0 

b 













the 


(which is 2 X 2 block tridiagonal) is a stationary G-chain in which p({l, 2}) = p({?7,, n-|-l}), n = 
2, 3, ■ ■ ■ is entangled. In this example p({n, n + j}), j > 2 is a product state with covariance 
matrix 

A/2 0 

0 A/2 


6. Conclusion 


An exchangeable chain of mean zero dnite mode Gaussian states is completely determined 
by two matrices A, B such that A — B is the covariance matrix of a Gaussian state and B is a 
nonnegative dednite matrix. Its asymptotic entropy rate is equal to the von Neumann 

entropy of the mean zero Gaussian state with covariance matrix A — B. 
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A class of examples of stationary Gaussian chains is constructed and their asymptotic entropy 
rates are evaluated. An example of a stationary entangled chain is presented. 
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